Raman scattering in Ceo and C48N12 aza-fullerene: First-principles study 
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We carry out large scale ab initio calculations of Raman scattering activities and Raman-active 
frequencies (RAFs) in C48N12 aza-fullerene. The results are compared with those of Ceo- Twenty- 
nine non-degenerate polarized and 29 doubly-degenerate unpolarized RAFs are predicted for C48N12. 
The RAF of the strongest Raman signal in the low- and high-frequency regions and the lowest and 
highest RAFs for C48N12 are almost the same as those of Ceo- The study of Ceo reveals the 
importance of electron correlations and the choice of basis sets in the ab initio calculations. Our 
best calculated results for C60 with the B3LYP hybrid density functional theory are in excellent 
agreement with experiment and demonstrate the desirable efficiency and accuracy of this theory for 
obtaining quantitative information on the vibrational properties of these molecules. 

PACS (numbers): 36.20.Ng, 31.15.Ar, 78.30.-j, 78.30.Na, 82.80.Gk 



I. INTRODUCTION 

In 1985, Ceo, a fascinating molecule formed as a trun- 
cated icosahedron with 20 hexagonal and 12 pentagonal 
faces, and 60 vertices, was discovered by Kroto et ed. [1]. 
Since then, a third form of pure carbon, called fullerenes 
[2], has been extensively studied. This kind of molecule 
can crystallize in a variety of three-dimensional struc- 
tures [3,4], with an even number of three-coordinated sp 2 
carbon atoms that arrange themselves into 12 pentagonal 
faces and any number (> 1) of hexagonal faces [2]. The 
macroscopic synthesis of soot [3] , which contains Ceo and 
other fullerenes in large compounds, plus the straightfor- 
ward purification techniques, which make pure fullerene 
materials available, opened new research opportunities in 
science, engineering and technology [5-9]. In the mean- 
time, doped fullerenes have also attracted a great deal of 
researchers' interest due to their remarkable structural, 
electronic, optical and magnetic properties [5-9]. For 
example, the doped tubular fullerenes can exhibit large 
third-order optical nonlinearities [8] and be ideal can- 
didates as photonic devices [7] for all-optical switching, 
data processing, eye and sensor protection (optical lim- 
iter). Another example is alkali metal-doped Cqq crys- 
tals, which can be superconducting [5] . Because of their 
unique structure and electronic properties [5], fullerenes 
can be doped in several different ways including endo- 
hedral doping [5,10], substitutional doping [11-16], and 
exohedral doping [5-7,9]. 

Since the average carbon-carbon bond length in Cgo is 
slightly larger than that in graphite, which can only be 
doped by boron, and the force constants [5] are some- 
what weakened by the curvature of the Cgo surface, both 
boron and nitrogen can substitute for one or more carbon 
atoms in Ceo [11-16]. In 1991, Smalley and coworkers 
[11] successfully synthesized boron-substituted fullerenes 
C6o-«B„ with n between 1 and at least 6. In 1995, Hum- 



melen et ed. [12] reported a very efficient method of syn- 
thesizing C59N, which has led to a number of detailed 
studies of the physical and chemical properties of CsgN 
[6,7,9,13]. Recently, Hultman et ed. [14] have synthesized 
nitrogen-substituted derivatives of Ceo with more than 
one nitrogen atom and reported the existence of a novel 
C48N12 aza-fullerene [14,15]. Very recently, we have stud- 
ied the bonding, electronic structure, Mulliken charge, in- 
frared (IR) spectrum, and NMR of C48N12 by using den- 
sity functional theory (DFT) and the 6-31G basis set [16]. 
We characterized 58 IR spectral lines, eight 13 C and two 
15 N NMR spectral signals of C48N12, and demonstrated 
that this aza-fullerene has potential applications as semi- 
conductor components for nanometer electronics because 
of its small energy gap, as a promising electron donor for 
molecular electronics, and as a good diamagnetic mate- 
rial because of the enhancement of diamagnetic factors in 
the carbon atoms [16]. The characterization of C48N12 is 
a timely problem both from the viewpoint of its practical 
applications and to understand doping in Cgo derivatives. 
In the present paper, we characterize the Raman spec- 
trum of C48N12. 

As a material is doped with foreign atoms, its mechan- 
ical, electronic, magnetic and optical properties change 
[5,7]. The ability to control such induced changes is vital 
to progress in material science. Raman and IR spectro- 
scopic techniques [17,18] are basic, useful experimental 
tools to investigate how doping modifies the structural 
and dynamical properties of the pristine material and to 
understand the physical origin of such induced changes. 
Over the past 10 years, both techniques have been widely 
used to study the vibrational properties of Ceo [19-24] , its 
derivative compounds [25,26], and (doped) carbon nan- 
otubes [27]. It has been shown that Ceo has in total 46 
vibrational modes including 4 IR-active [19-22] and 10 
Raman-active [23,24] vibrational modes. Well-resolved 
Raman spectra [25] are also available for Ceo and a num- 
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ber of its derivative compounds. 

In this paper, we perform first-principles calculations 
of Raman scattering activities (RSAs) and Raman-active 
frequencies (RAFs) in both C48N12 and Cgo using density 
functional theory (DFT) and the restricted Hartree-Fock 
(RHF) method. Very recently Choi et al. [28] have theo- 
retically assigned all 46 vibrational modes of Cgo , includ- 
ing a scaling of the force field by using Pulay's method. 
In this paper, however, we carry out ab initio a series of 
calculations for Cw We want to test the efficiency and 
accuracy of such first-principles calculations and study, 
in detail, basis set effects on RSAs and RAFs by com- 
paring our theoretical results with available experiments 
[5,23-26]. To the best of our knowledge, such basis set ef- 
fects on the RSAs and RAFs of Ceo have not been consid- 
ered before. These calculations for Ceo give us a bench- 
mark for assigning the Raman-active vibrational modes 
of C48N12 and provide us with constructive insight into 
the microscopic mechanisms responsible for the difference 
between C48N12 and Ceo- We find that the 10 RAFs for 
Ceo, obtained by using a hybrid DFT method and large 
basis sets, are in excellent agreement with Raman ex- 
periments. We predict that C48N12 aza-fullerene has 58 
RAFs including 29 non-degenerate polarized modes and 
29 doubly-degenerate unpolarized modes. 

This paper is organized as follows. Section II briefly 
reviews the ab initio methods, basis sets and the the- 
ory of calculating Raman scattering activities. Section 
III presents our Raman results for both Ceo and C48N12 
obtained by using RHF and DFT methods, and the re- 
sults are compared to results obtained by other theoret- 
ical methods. Our conclusions are given in section IV. 



II. THEORY 

A. Ab Initio Methods and Basis Sets 

Ab initio methods obtain information by solving vari- 
ationally the Schrodingcr equation without fitting pa- 
rameters to experimental data. Instead, experimen- 
tal data guides the selection of the ab initio methods 
rather than directly entering the computational proce- 
dures. Fullerenes have been challenging molecules for 
ab initio calculations because of their size [29]. Recent 
advances in ab initio electronic structure methods and 
parallel computing have brought a substantial improve- 
ment in the capabilities to predict and study the proper- 
ties of large molecules. The coupled cluster method [30] 
has been used to predict phenomena in C 2 o [31]. Other 
ab initio methods, which are less demanding in terms of 
computation cost than the coupled cluster method, have 
been used for much larger fullerenes and carbon nan- 
otubes, for example, Ceo [32,33] with self-consistent field 
and Moller-Plcssct second-order (MP2) theory, C240 [34] 
and carbon nanotubes [35] with density functional the- 
ory (DFT), and C540 [36] with the Hartree-Fock (HF) 



method. The major expense in HF and DFT calcula- 
tions arises from solving the electronic quantum Coulomb 
problem. The effective Hamiltonian diagonalization (a 
procedure that scales as Njf, Nb being the number of 
basis functions) represents only a minor portion of the 
computational time in calculations of molecular clusters 
containing up to several hundred atoms [37]. In addi- 
tion, DFT [38] requires an additional three-dimensional 
numerical quadrature to obtain the exchange and corre- 
lation energies [30]. Both HF and DFT methods have 
been implemented into the Gaussian 98 program [39] . In 
this paper, all calculations are performed by using the 
Gaussian 98 program [39,40]. 

One of the approaches inherent in all ab initio methods 
is the introduction of a basis set [30] . If the basis set is 
complete, exact expansions of the molecular orbitals can 
be obtained. However, a complete basis set requires an 
infinite number of functions, which is impossible in ac- 
tual calculations. Generally, a smaller basis set provides 
a poorer representation. Moreover, only the parts of the 
molecular orbital which correspond to the selected basis 
can be represented. Since the computational effort of ab 
initio methods scales formally as N£ [30], it is impor- 
tant to make the basis set as small as possible without 
compromising the accuracy or missing part of the state 
space which should be represented. Hence, one purpose 
of this paper is to determine the effects of basis sets [39] 
for the calculations of Raman spectra by consdering the 
Slater-type basis set STO-3G and split valence basis sets 
3-21G, 6-31G and 6-31G(d) (i.e., 6-31G*). 



B. Vibrational Analysis and Raman Scattering 
Activity 

The vibrational analysis of polyatoms described by 
Wilson et al. [41] has been implemented in the Gaussian 
98 program. This analysis is valid only when the first 
derivatives of the energy with respect to the displacement 
of the atoms are zero (in other words, the geometry used 
for vibrational analysis must be optimized at the same 
level of theory and with the same basis set that the sec- 
ond derivatives are generated with). The force constant 
matrix K is defined as the second partial derivatives of 
the potential V with respect to the displacement of the 
atoms in cartesian coordinates (for example, Axk, Ay^, 
Azk of the fcth atom), i.e., a in x 3n matrix (n is the 
number of atoms) whose elements are given by 

where (771, 772, 773) = (Aari, Aj/i, A21). (...)o means that 
the second partial derivatives are taken at the equilib- 
rium positions of the atoms. Usually, the matrix K in 
cartesian coordinates is converted to a new matrix K in 
mass-weighted cartesian coordinates % = m^ 2 ?,,, i.e., 
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where mi is the mass of the atom. Then, the eigenval- 
ues Xk of K give the fundamental frequencies (in the 

unit of cm -1 ), i.e., Vk = AjJ/ 2 / (2irc), where c is the veloc- 
ity of light in vacuum. The eigenvectors give the normal 
modes. 

To lowest order, Raman intensities are proportional to 
the derivatives of the dipole polarizability with respect 
to the vibrational normal modes of the material, eval- 
uated at the equilibrium geometry. For example, most 
Raman scattering experiments use a plane-polarized in- 
cident laser beam. The direction of the incident beam, 
the polarization direction of this beam and the direction 
of observation are chosen to be perpendicular to each 
other. Under these circumstances, the first-order differ- 
ential cross section for Raman scattering in the gth vi- 
brational mode is written as [42] 
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Iraman is the Raman scattering activity, and n q = 
cxp [— TiiOq/ (kbT)]. lo s is the frequency of the scattered 
radiation, ui q is the frequency of the gth vibrational mode, 
«b is Boltzmann's constant, and T is temperature. For 
the special case of ir/2 scattering geometry, I ra man can be 
written as 
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where S g is the normal coordinate corresponding to the 
qth vibrational mode and a is the dipole polarizability 
tensor. 

To obtain the Raman activities, one must compute the 
derivatives of the polarizability with respect to the nor- 
mal mode coordinates S g . These can be viewed as direc- 
tional derivatives in the space of 3N nuclear coordinates 
and are expressed in terms of derivatives with respect to 
atomic coordinates, Rk- For the polarizability compo- 
nent aij (i, j = x, y, z), we have 
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where £kq = dRk/dEq is the fcth atomic displacement of 
the qth normal mode. Then, the necessary derivatives 



can be expressed in terms of the atomic forces as follows 
[41,42] 
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where E is the total energy, Gi is the «th component of 
an assumed external electric field G, and F k = dE/dRk 
is the calculated force on the fcth atomic coordinate. 



III. RESULTS 

A. Geometry Optimizations 

In our previous work [16], we calculated the bond 
lengths of C 48 Ni2 by using the B3LYP [43] hybrid DFT 
method with the 6-31G basis set. The geometry of 
C48N12 is shown in Fig.l, where the 10 unique sites (ver- 
tices 1 to 10) can be identified from nuclear magnetic 
resonance [16]. It was found that there are 15 unique 
types of bonds in C48N12 [16]: 6 nitrogen-carbon bonds 
and 9 carbon-carbon bonds. In contrast, Cqq has one 
kind of single (C-C) bond and one kind of double bond 
(C=C). 




FIG. 1.: C48N12 structure. Red is for nitro- 
gen sites and grey for carbon sites. There are 
10 unique vertices (1 to 10) as labeled. The 15 
bonds between labeled vertices are all unique. 

In this paper, we discuss the effect of basis sets using 
different ab initio methods. Before calculating the RSAs 
and RAFs of both Cgo and C48N12, we first optimize 
their geometries by using RHF and B3LYP hybrid DFT 
methods with a variety of basis sets including STO-3G, 
3-21G, 6-31G and 6-31G*. The B3LYP DFT method in- 
cludes a mixture of Hartree-Fock (exact) exchange, Slater 
local exchange [44], Becke 88 non-local exchange [45], 
the VWN III local exchange-correlation functional [46] 
and the LYP correlation functional [47]. The optimized 
bond lengths for C48N12 and Ceo are listed in Table I. 
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Compared to Ceo, in C48N12 the lengths of bonds that 
were Ceo double bonds (for example, C1-C2 , C3-C10) in- 
crease, while the lengths of bonds that were single bonds 
in Ceo (for example, C3-C4, C6-C10) decrease, yielding 
two bonds with a small difference in length. According 
to our B3LYP/6-31G* calculation, the bond length dif- 
ference is reduced from 0.059 A in C 60 to 0.027-0.041 A 
in C48N12. This can be explained by the fact that the 12 
extra electrons in C48N12 fill t\ u and t\ g orbitals which 
have some antibonding character and consequently in- 
crease the length of C=C double bonds. Compared with 
B3LYP calculations and experiments [48] for Ceo, RHF 
generally underestimates the lengths of C=C and C-C 
bonds by about 0.5%. Judging from the success of the 
B3LYP geometries for Ceo, we conclude that the role 
of electron correlation is important in the accurate de- 
scription of localization/delocalization of 7r electrons in 
Ceo- Such bond equalization has also been observed in 



the DFT calculations for K3C60 by Bohnen ct al. [49] and 
for C^q 6 by Choi et al. [28] . To prevent the collapse of the 
valence functions into the inner shell and ensure an ad- 
equate description of bonding interactions which involve 
overlap of valence funcitons, the 6-31G* basis set for car- 
bon and nitrogen provides a better description of the 
inner-shell region as well as the valence region than STO- 
3G, 3-21G and 6-31G basis sets [39,62]. Overall, the bond 
length L in Ceo or C48N12, as listed in Table I, follows the 
order: L 6 _ 3 ig* < L 3 _ 2 ig < L 6 _ 3 i G < L S to-3G, where 
L6-31G* for C 6 o is in good agreement with experiment. 
The ordering of bond lengths suggest that increasing the 
flexibility of the basis from 3-21 to 6-3 1G favors expan- 
sion of the Ceo, while the extra polarization functions in 
6-31G* are needed to recontract Ceo- As a result, the 
3-2 1G basis set fortuitously gives bond lengths closer to 
experiment than does the bigger 6-31G basis set. 



Table I: Bond lengths (L, in A ) in C48N12 and Cqq calculated by using B3LYP hybrid DFT and RHF methods 
with STO-3G, 3-21G, 6-31G and 6-31G* basis sets, i and j in Ci and Nj are the site numbers labeled in Fig.l. The 
symbols "C-C" and "C=C" in parenthesis denote the original single and double bonds in Cqo- The averaged bond 
length L ave = (L c =c + 2L C -c)/3. 



STO-3G 3-21G 6-31G 6-31G* EXP. 
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Ref. [48] 
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.417 
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.408 


1.422 
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.413 


1.416 


1.402 
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.428 
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.391 
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.386 


1.413 


1 


.388 


1.406 


1.384 
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.474 
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.460 
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.432 
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.421 


1.432 
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.416 


1.430 


1.427 
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.470 
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.466 
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.452 
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.452 


1.449 


1 


.448 


1.446 


1.452 




1 


.409 


1 


.359 


1 


.395 
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.362 


1.397 


1 


.365 


1.390 


1.359 
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.456 
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.457 
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.432 
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.435 


1.434(6) 
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.437 


1.431(4) 


1.436(8) 
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.457 
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.433 
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.437 


1.435(4) 
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1.437(2) 
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.455 
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.431 
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.413 


1.429 
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.411 


1.422 
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.464 


1 


.430 
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.429 
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1.429 


1.431 
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1.394 
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1.415 
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C 48 Ni2 C1-C2 (C=C) 
C1-C2 (C-C) 
C1-N5 (C-C) 
C3-C2 (C-C) 
C3-C4 (C-C) 
C3-C10 (C=C) 
C4-C8' (C=C) 
C4-N5 (C-C) 
N5-C6' (C=C) 
C6-C10 (C-C) 
C6-C7 (C-C) 
C7-C8 (C-C) 
C7-N9' (C=C) 
C8-N9 (C-C) 
N9-C10 (C-C) 
C 60 C=C (C1-C2') 



Lave 



According to our DFT or RHF calculations, the di- 
ameter of C48N12 varies, as compared to that of Ceo, 
from -0.3 A to +0.2 A, which has a strong effect on the 
force constant that determines the frequencies of the vi- 
brational modes of molecule. Hence, we expect that the 
vibrational frequencies of C48N12 will be strongly influ- 
enced by the bond equalization effects and the expansion 
and contraction of the ball. In the following, we investi- 
gate the vibrational properties of C48N12 and Ceo based 
on the optimized geometries. 



B. Raman-active Vibrational Frequencies 

Using the Gausian 98 program [39,40], we calculated 
the RAFs of both Ceo and C48N12. We compare with 
Ceo experiment to determine the effects of electron cor- 
relations and basis sets. It should be mentioned that our 
frequencies have not been scaled. 

Table II summarizes the calculated results for Ceo ob- 
tained by using RHF and B3LYP hybrid DFT methods. 
As shown by Dresselhaus et al. [50], there are 46 dif- 
ferent vibrational modes in the 174 independent normal 
vibrations of Ceo • These modes are classified in even and 
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odd parities and in the ten irreducible representatives 
of the I h point group [50]: the {a g , a u }, {t lg , t\ u , t 2gi 
t2u}, {9g, 9u\ and {h g , h u } modes are non-, threefold- 
, fourfold- and fivefold-degenerate, respectively. Among 
the 46 vibrational modes are 10 Raman-active ones. In 
choosing a basis set for the first-principles calculation, 
one must make a compromise between accuracy and CPU 
time. Without significant computational cost, one can 
do B3LYP/STO-3G calculation and still obtain results 
more accurate than any RHF calculations. Going be- 
yond STO-3G for B3LYP cases requires a drastic increase 
in CPU time. Suprisingly, going just to 3-21G provides 
the most accurate results, while for the bigger basis set 
6-31G, the results are worse and adding a polarized func- 
tion to 6-31G only slightly improves the results. The 3- 
21G basis set gives systematically lower frequencies than 
the 6-3 1G basis set, while the frequencies obtained from 



the 6-31G* basis set typically lie in between the results 
of the other two basis sets. In contrast, 6-31G* does pro- 
vide the most accurate bond lengths. This suggests that 
the better accuracy of 3-21G is foirtuitous. Increasing the 
basis set to 6-31G stiffens the bonds, while adding the po- 
larization function compensates by softening the bonds. 
In comparison with the B3LYP results, RHF calculated 
frequencies are too high due to an incorrect description of 
bond dissociation which leads to an increased force con- 
stant, while B3LYP's with large basis sets are generally 
in good agreement with the experiments of Wang et al. 
[21]. In comparison with both effects, we find that the 
correlation effect on the RAFs is stronger than the basis 
set effect. This demonstrates the importance of electron 
correlation in the accurate description of the vibrational 
frequencies of these molecules. 



Table II: RHF and B3LYP hybrid DFT calculations of Raman scattering activities ( I ra man, hi 10~ 14 m 4 /kg ) 
of Ceo with the corresponding vibrational modes and frequencies v (cm -1 ). h g and a g modes arc unpolarized and 
polarized, respectively. Numbers in the parenthesis are the relative errors between the calculated and the experimental 
frequencies v exp (see Table III) of Wang et al. [23] . 



Method Mode 




STO-3G 




3-21G 




6-31G 




6-31G* 


Iranian. 


V 


Iranian 


V 


Iranian 


V 


Iranian 


V 


RHF a g 


3894 


1684 (14.6%) 


3467 


1604 (9.2%) 


3840 


1637 (11.4%) 


3968 


1600 (8.9%) 




441 


553 (21.1%) 


488 


518 (5.0%) 


558 


526 (6.8%) 


574 


527 (6.4%) 




255 


1912 (21.6%) 


224 


1772 (12.7%) 


279 


1799 (14.4%) 


283 


1791 (13.9%) 




34 


1658 (16.3%) 


35 


1546 (8.4%) 


38 


1585 (11.1%) 


22 


1562 (9.5%) 




188 


1482 (18.7%) 


215 


1326 (6.2%) 


199 


1377 (10.3%) 


210 


1380 (10.6%) 




87 


1290 (17.4%) 


76 


1184 (7.7%) 


44 


1208 (9.9%) 


105 


1208 (9.9%) 




38 


886 (14.6%) 


38 


828 (7.2%) 


5 


843 (9.0%) 


46 


840 (8.7%) 




11 


836 (18.1%) 


16 


761 (7.5%) 


5 


821 (15.9%) 


9 


794 (12.2%) 




14 


509 (18.0%) 


19 


475 (10.2%) 


16 


496 (14.9%) 


12 


482 (11.9%) 




75 


302 (11.9%) 


111 


295 (9.2%) 


116 


296 (9.5%) 


115 


289 (7.1%) 


B3LYP a g 


3008 


1549 (5.4%) 


2643 


1501 (2.1%) 


2758 


1524 (3.7%) 


2760 


1504 (2.4%) 




654 


502 (1.9%) 


681 


491 (0.4%) 


750 


496 (0.6%) 


724 


489 (0.8%) 


h g 


320 


1677 (6.6%) 


297 


1609 (2.3%) 


325 


1627 (3.4%) 


312 


1618 (2.8%) 




28 


1500 (5.2%) 


18 


1436 (0.7%) 


23 


1466 (2.8%) 


19 


1455 (2.0%) 




138 


1332 (6.7%) 


168 


1231 (1.4%) 


151 


1274 (2.1%) 


170 


1275 (2.1%) 




74 


1166 (6.1%) 


71 


1112 (1.2%) 


75 


1129 (2.8%) 


71 


1125 (2.4%) 




47 


802 (3.8%) 


48 


781 (1.1%) 


44 


788 (1.9%) 


48 


766 (0.9%) 




4 


734 (3.7%) 


5 


678 (4.3%) 


4 


738 (4.3%) 


5 


718 (1.4%) 




4 


449 (4.2%) 


5 


429 (0.5%) 


4 


448 (3.9%) 


5 


436 (1.2%) 




95 


271 (0.5%) 


128 


271 (0.3%) 


135 


272 (0.2%) 


130 


266 (1.7%) 



For comparison, Table III lists the calculated vibra- 
tional frequencies of Cqo obtained by using various the- 
ories, for example, the semi-empirical MNDO [51] and 
QCFF/PI [52] methods. Of these, the QCFF/PI method, 
which has been parameterized mainly with respect to 
vibrational frequencies of conjugated and aromatic hy- 
drocarbons [53], results in the best results although it 
gives less satisfactory geometry. Such accurate predic- 
tion implies that the electronic structures of Cqq is not 
much different from other aromatic hydrocarbons [52]. 



Haser et al. [33] showed that the approximate harmonic 
frequencies for the two a g vibrational modes of Cgo are 
1615 cm -1 (9.9%) and 487 cm" 1 (1.2%) at HF/DZP, 
1614 cm- 1 (9.9%) and 483 cm" 1 (2.0%) at HF/TZP, 1614 
(9.8%) cm" 1 and 437 (12.2%) cm" 1 at MP2/DZP, and 
1586 (7.9%) cm" 1 and 437 (12.2%) cm" 1 at MP2/TZP, 
where the percentages in the parenthesis are the rela- 
tive errors of the calculated results to the experimental 
frequencies obtained by Wang et al. [23]. Their HF cal- 
culations are in agreement with our RHF/3-21G results. 
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Their MP2 results are more accurate when obtained with of electron correlation in predicting accurately the vibra- 
large basis sets, which also demonstrates the importance tional frequencies of a molecule. 

Table III: The vibrational frequencies obtained by other theoretical calculations performed by Choi et al. [28], 
Bohnen et al. [49], Stanton et al. [51], Negri et al. [52], Jishi et al. [57], Dixon et al. [58], Hara et al. [59] and Onida et 
al. [60]. Numbers in the parenthesis are the relative errors to the experimental frequencies, v exp , of Wang et al. [23]. 



SIFC 


LDA 


LDA 


LDA 


CPMD 


QCFF/PI 


MFCM 


MNDO 


Exp. 


[28] 


[59] 


[58] 


[49] 


[60] 


Ref. [52] 


[57] 


[51] 


[23] 



Clg 

1474 (0.3%) 1531 (4.2%) 1525 (3.8%) 1475 (0.4%) 1447 (1.5%) 1442 (1.8%) 1468 (0.1%) 1667 (13.5%) 1469 

484 (1.8%) 502 (1.8%) 499 (1.2%) 481 (2.4%) 482 (2.2%) 513 (4.1%) 492 (0.2%) 610 (23.7%) 493 

h g 

1582 (0.6%) 1609 (2.3%) 1618 (2.8%) 1580 (0.4%) 1573 (0.0%) 1644 (4.5%) 1575 (0.1%) 1722 (9.5%) 1573 

1419 (0.5%) 1475 (3.4%) 1475 (3.4%) 1422 (0.3%) 1394 (2.2%) 1465 (2.7%) 1401 (1.8%) 1596 (11.9%) 1426 

1250 (0.2%) 1288 (3.2%) 1297 (3.9%) 1198 (4.0%) 1208 (3.2%) 1265 (1.4%) 1217 (2.5%) 1407 (12.7%) 1248 

1117 (1.6%) 1129 (2.7%) 1128 (2.6%) 1079 (1.8%) 1098 (0.1%) 1154 (5.0%) 1102 (0.3%) 1261 (14.7%) 1099 

782 (1.2%) 794 (2.7%) 788 (1.9%) 763 (1.3%) 775 (0.3%) 801 (3.6%) 788 (1.9%) 924 (19.5%) 773 

704 (0.6%) 711 (0.4%) 727 (2.7%) 716 (1.1%) 730 (3.1%) 691 (2.4%) 708 (0.0%) 771 (8.2%) 708 

436 (1.2%) 430 (0.2%) 431 (0.0%) 422 (2.1%) 435 (0.9%) 440 (2.1%) 439 (1.9%) 447 (3.7%) 431 

272 (0.7%) 269 (0.4%) 261 (3.3%) 263 (2.6%) 261 (3.3%) 258 (4.4%) 269 (0.4%) 263 (2.7%) 270 



In addition, a number of force-constant models (FCMs) 
that include the interactions up to the second-nearest 
neighbors [54-56] are used to calculate the vibrational 
frequencies of Cqq . None of them yield good agreement 
with the experiment data. For example, an empirical 
force field, which has been parameterized with respect 
to polycyclic aromatic hydrocarbons, is used with the 
Hiickcl theory and predicts that the vibrational frequen- 
cies of the two a g modes are 1409 cm' 1 and 388 cm -1 
[55], which are too low. However, the modified FCM 
(MFCM) by Jishi et al. [57] considered interactions up 
to the third-nearest neighbors. The results obtained by 
using MFCM [57], as shown in Table III, are in excellent 
agreement with the experiments of Wang et al. [21]. 

Table III also list the calculated vibrational frequen- 
cies of Cqo obtained by other DFT methods, for example, 
local density approximation (LDA) [49,58,59] and DFT- 
LDA-based Car-Parrinello molecular dynamics (CPMD) 
simulation [60]. In general, those calculated results are 
in good agreement with experiment. 

Recently, Choi et al. [28] have performed B3LYP vi- 
brational calculations of Cgo with a 3-2 1G basis set but 
involving scaling of the internal force constants (SIFC) 
K™* by using Pulay's method [61], i.e., 

^scaled = ( a . Sj )l/2£jrrt (7) 

where K™* is the force constant in internal coordinates 
( the Gaussian 98 program [39] uses this form), and 
and Sj are scaling factors for the zth and jth redundant 
internal coordinates, respectively. They optimized the 
scaling factors by minimizing the root-mean-square de- 
viations between the experimental and calculated scaled 
frequencies. Their results are also listed in Table III. 



Their scaling procedure does improve the accuracy for 
the 10 Raman-active vibrational frequencies of Cgo- 

It should be mentioned that although vibrational cal- 
culations are all done for an isolated Cqo molecule, the 
calculated results are usually compared to experimental 
data on Ceo molecular vibrations in the solid state, in the 
gas-phase or in solution. The theoretical justification for 
this lies in the weak Van der Waals bonds between 
molecules [5] . The comparison between the observed vi- 
brational spectra of C@o [5] demonstrated that the mea- 
sured mode frequencies are slightly different from one 
phase to another, showing that the interaction between 
Ceo molecules is weak. 

In Table IV and V, we list the RAFs for C 48 Ni 2 cal- 
culated with RHF and B3LYP hybrid DFT methods , 
respectively, and different basis sets. In contrast with 
Ceo, we find that there are in total 116 different vibra- 
tional modes [16] for C48N12 because of its lower symme- 
try, Sq [16]. These vibrational modes are classified in 58 
doubly-degenerate and 58 nondegenerate modes. Among 
those vibrational modes, there are 58 IR-active vibra- 
tional modes [16] and 58 Raman-active modes includ- 
ing 29 doubly-degenerate and 29 non-degenerate ones as 
listed in Table IV and V. Similar to Ceo, including the 
electron correlation and increasing the basis size leads to 
a redshift of the RAFs of C48N12. From our calculated 
results, we find that the nitrogen-substitutional doping 
results in a symmetry lowering and an increase in the re- 
duced mass. The symmetry lowering splitts some of the 
degenerate vibrational modes observed in Ceo and makes 
many more modes Raman-active for C48N12. Overall, 
the increase of the reduced mass red-shifts the vibrational 
frequencies of C50 . 
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C. Raman Scattering Activities 

In this paper, we only calculate non-resonant Raman 
intensities. We performed ab initio calculations of Ra- 
man scattering activities I r aman [17,18] for the optimized 
geometries of C48N12 and Ceo by using the Gaussian 
98 program [39,40] with RHF and B3LYP hybrid DFT 
methods. 



Doubly-degenerate Modes 



STO-3G 3-21G 6-31G 6-31G* 



Table IV: Fifty eight Raman-active frequencies v 
(cm _1 )for C48N12 calculated by using RHF method with 
STO-3G, 3-21G, 6-31G and 6-31G* basis sets. 



Doubly-degenerate Modes 



STO-3G 3-21G 6-31G 6-31G H 



262.9 
285.8 
429.1 
470.4 
513.6 
583.2 
603.5 
651.8 
725.4 
773.4 
780.2 
831.2 
848.0 
870.7 
934.6 
961.3 
1238.4 
1280.9 
1313.3 
1422.9 
1476.0 
1546.2 
1608.0 
1647.0 
1716.5 
1780.9 
1831.0 
1882.8 
1915.5 



262.1 

280.9 

408.9 

433.3 

478.7 

516.9 

609.5 

647.3 

698.9 

726.4 

761.6 

775.8 

820.3 

843.6 

942.6 

974.4 

1159.9 

1196.7 

1238.4 

1276.3 

1340.6 

1420.7 

1468.3 

1490.0 

1541.4 

1631.3 

1669.5 

1707.7 

1735.1 



268.3 

285.5 

427.5 

454.6 

496.7 

548.6 

611.3 

649.5 

736.4 

759.2 

785.1 

806.9 

833.4 

852.5 

943.5 

965.8 

1191.9 

1225.0 

1261.7 

1324.6 

1382.7 

1456.9 

1517.1 

1537.2 

1591.4 

1670.8 

1709.9 

1749.3 

1777.7 



260.7 

280.8 

412.9 

440.1 

488.3 

539.9 

592.2 

628.6 

706.2 

740.9 

757.4 

788.0 

820.8 

840.1 

917.9 

932.8 

1164.4 

1214.5 

1249.9 

1317.6 

1377.4 

1448.6 

1503.4 

1540.1 

1586.5 

1663.5 

1719.9 

1758.7 

1788.1 



Non-degenerate Modes 



STO-3G 3-21G 6-31G 6-31G* 



293.9 
444.6 
462.5 
510.2 
522.6 
575.6 
608.3 
645.7 
650.4 
723.1 
744.2 
777.3 
831.6 
866.8 
923.6 
952.9 
1216.2 
1327.5 
1344.1 
1400.2 
1468.7 
1498.2 
1533.9 
1583.4 
1655.1 
1721.2 
1766.4 
1817.4 
1946.3 



287.6 

406.1 

415.1 

455.2 

490.8 

520.8 

566.1 

637.1 

650.6 

654.8 

705.2 

759.7 

778.6 

824.8 

918.4 

964.6 

1164.9 

1243.7 

1259.7 

1335.0 

1347.1 

1374.6 

1414.9 

1436.5 

1527.5 

1559.5 

1612.4 

1643.5 

1766.7 



291.2 

428.9 

447.6 

483.2 

511.4 

542.8 

586.8 

639.7 

652.2 

670.3 

748.8 

783.6 

802.4 

836.1 

919.0 

955.7 

1192.3 

1264.1 

1281.3 

1374.0 

1389.0 

1424.4 

1468.7 

1485.1 

1578.0 

1612.3 

1648.6 

1686.6 

1803.8 



287.9 
387.5 
424.5 
454.0 
488.6 
523.7 
575.2 
614.7 
636.8 
651.8 
716.4 
738.3 
772.8 
819.0 
897.0 
916.2 
1174.7 
1251.1 
1272.8 
1347.6 
1376.6 
1380.0 
1429.1 
1475.9 
1551.9 
1602.1 
1625.7 
1680.1 
1808.7 



245.1 

258.9 

376.1 

410.2 

446.7 

492.6 

548.8 

580.3 

617.3 

649.4 

665.2 

715.9 

766.6 

781.7 

840.8 

858.1 

1123.7 

1162.2 

1181.5 

1274.1 

1318.1 

1386.2 

1451.9 

1466.1 

1503.9 

1570.8 

1595.2 

1633.1 

1661.7 



247.7 

261.1 

367.7 

388.2 

428.7 

449.5 

567.5 

596.3 

612.8 

629.1 

670.6 

695.0 

760.4 

776.0 

865.3 

887.4 

1084.4 

1122.8 

1146.8 

1186.0 

1243.4 

1310.8 

1369.3 

1387.9 

1415.3 

1488.6 

1516.0 

1544.0 

1577.6 



252.1 

264.5 

386.7 

408.9 

446.2 

482.0 

566.0 

596.8 

652.2 

662.4 

698.1 

713.9 

768.0 

782.4 

860.2 

876.6 

1104.9 

1140.6 

1162.9 

1220.0 

1273.6 

1336.6 

1406.8 

1420.4 

1451.6 

1516.8 

1550.6 

1577.6 

1609.8 



245.1 

259.5 

371.2 

396.0 

436.5 

471.6 

551.1 

580.6 

626.6 

641.2 

671.1 

698.6 

765.8 

779.7 

842.6 

853.9 

1092.7 

1138.0 

1155.8 

1222.6 

1271.7 

1336.3 

1404.1 

1419.7 

1451.0 

1515.4 

1551.0 

1577.8 

1609.8 



Non-degenerate Modes 



STO-3G 3-21G 6-31G 6-31G* 



263.8 

381.9 

405.7 

444.4 

471.2 

498.3 

544.1 

577.5 

585.2 

621.2 

624.9 

657.8 

729.5 

768.2 

835.5 

855.9 

1114.7 

1193.6 

1220.7 

1315.9 

1337.1 

1366.7 

1416.4 

1438.7 

1481.0 

1530.7 

1558.4 

1580.5 

1682.7 



264.3 
368.2 
397.7 
415.4 
458.4 
490.5 
504.5 
588.2 
597.8 
603.0 
614.1 
657.1 
690.1 
764.8 
850.5 
881.0 
1079.8 
1160.0 
1181.4 
1240.9 
1253.2 
1292.5 
1346.6 
1372.8 
1418.8 
1440.4 
1487.1 
1508.5 
1592.0 



267.5 

389.4 

415.0 

441.7 

477.3 

500.2 

544.2 

592.1 

601.5 

616.3 

652.1 

684.2 

709.4 

772.5 

844.4 

868.6 

1099.8 

1172.5 

1195.8 

1267.3 

1280.2 

1322.1 

1382.6 

1414.5 

1453.2 

1477.2 

1511.4 

1545.1 

1623.7 



263.6 

376.2 

398.3 

423.8 

466.6 

494.9 

510.0 

575.7 

588.1 

596.9 

627.0 

656.2 

688.5 

765.5 

829.7 

843.7 

1092.1 

1161.6 

1189.0 

1265.1 

1279.0 

1319.8 

1378.9 

1397.4 

1441.1 

1476.6 

1505.3 

1530.4 

1622.9 



Table V: Fifty eight Raman-active frequencies v (cm 1 )for 
C48N12 calculated by using B3LYP hybrid DFT methods with 
STO-3G, 3-21G, 6-31G and 6-31G* basis sets. 



To test the basis set dependence of the theoretical 
RSAs, calculations for basis sets of different sizes have 
been performed for Ceo- In Table II are listed the RSAs 
for C 60 calculated with STO-3G, 3-21G, 6-31G and 6- 
31G*. For each basis set, the strongest Raman-active 
lines of Ceo are the two a g modes, which is consistent with 
Raman experiments [20-24] . The a g modes are identified 
by their polarized character, which strongly suggests that 
both a g modes are totally symmetric [50]. The remain- 
ing 8 Raman-active modes are unpolarized, consistent 
with the fivefold-degenerate h g symmetry [50]. Also, we 
find that electron correlation and basis set change signifi- 
cantly the RSA of a given Raman-active mode. However, 
the correlation effect is stronger than the basis set effect. 
In general, electron correlation reduces the RSA of the 
strongest Raman-active a g mode of Ceo by about 25%. 
For other specific Raman-active modes, the electron cor- 
relation predicts at least a 20% decrease (or increase) of 
the RSA obtained by the RHF method. Increasing the 
size of the basis set from STO-3G to 6-31G results in 
about a 10% decrease (or increase) of RSA, while adding 
polarization functions leads to about a 3% decrease (or 
increase) of RSA. For a specific mode (for example, the 
gth one), we find that the fcth atomic displacement t;kq of 
this mode changes little as we choose different basis sets 
or/and consider the electron correlation. Thus, accord- 
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ing to Eq.(5), the RSA mainly depends on the derivative 
of the polarizability a with respect to atomic coordinates 
R, which are related to the choice of basis sets and the 
inclusion of electron correlations. 



0.1 



(a) C, 8 N 12 : polarized □ 

v ' C 48 N 12 : unpolarized • 

C 60 : polarized 
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FIG. 2: Ah initio calculations of Raman scat- 



tering activities (I ra r, 



in 10- 14 m 4 /kg) 



C 48 Ni2 with (a) RHF/3-21G and (b) B3LYP/3- 
21G. Open squares and filled circles are non- 
degenerate polarized and doubly-degenerate un- 
polarized Raman-active modes, respectively for 
C48N12. The solid and dot-dashed lines are 
the calculated unpolarized and polarized Raman 
spectral lines of Ceo, respectively. 




FIG. 3: The vibrational displacements of 
sites 1 to 5 for the strongest Raman spec- 
tral lines in both low- and high-frequency re- 
gions for the B3LYP/3-21G case. The open 
and filled circles are for the low- and high- 
frequency cases of Ceo, respectively. The 
open and filled triangles arc for the low- and 
high-frequency cases of C48N12, respectively. 

In Fig. 2 (a) (b), we show the calculated RSAs I ram an at 
the corresponding RAFs v for C48N12 by using RHF/3- 
21G and B3LYP/3-21G, respectively. We find that the 
Raman spectrum of C48N12 separates into two regions, 
i.e., high-frequency (1100 cm -1 to 1700 cm -1 ) and low- 
frequency (200 cm -1 to 1000 cm -1 ) regions, which are 
similar to those of Ceo- In detail, this aza-fullerene, un- 
like C 60 , has an equal number of polarized and unpolar- 
ized Raman-active modes in each region and, in particu- 
lar, has 6 more Raman-active modes in the low-frequency 
region than in the high-frequency one. The strongest 
Raman spectral lines in both low- and high-frequency re- 
gions are the non-degenerate polarized ones located at 
1509 cm -1 and 491 cm -1 , respectively, which are almost 
the same as those of the strongest a g modes of Ceo- Tak- 
ing B3LYP/3-21G calculations as an example, we show in 
Fig. 3 the vibrational displacements of sites 1 to 5 for the 
strongest Raman modes in both C48N12 and Ceo- For 
Ceo, the pentagon structure for the low- frequency case 
contracts largely in the x-y plane and shows strong col- 
lective z-axis vibrations, while the pentagon structure for 
the high-frequency case expands largely in the x-y plane 
and shows small collective z-axis vibrations. In contrast, 
as a result of N doping, the pentagon structure for C48N12 
for the low-frequency case expands slightly and shows col- 
lective vibration along the z-x direction, while the pen- 
tagon structure for the high-frequecny case contracts, ac- 
companying a large stretching of site 4. As shown in 
Fig. 2, the lowest and highest RAFs of C48N12 are al- 
most the same as those of Cgo- Comparing Fig. 2 (a) with 
Fig.2(b) shows that the RAFs of C 48 Ni2 predicted by 
RHF method, as discussed above, are redshifted due to 
electron correlations. The effects of basis set and electron 
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correlations on RSAs of C48N12 are similar to those ob- 
tained for Cgo • Detailed analysis of Raman-active vibra- 
tional modes of C48N12 shows that (i) the Raman spectra 
in the high-frequency region (> 1400cm -1 ) comes mainly 
from carbon-carbon vibrations; (ii) the strong nitrogen- 
carbon vibrations occur only in the low-frequency region; 
(iii) the strongest and weakest Raman signals come from 
the contributions of carbon-carbon vibrations. 



IV. SUMMARY 

In summary, we have performed large scale ab initio 
calculations of RSAs and RAFs in C48N12 as well as Ceo 
using B3LYP hybrid DFT and RHF methods. We pre- 
dict that C4 8 N 12 has 29 non-degenerate polarized and 
29 doubly-degenerate unpolarized RAFs, and the RAF 
of the strongest spectral lines in the low- and high- fre- 
quency regions and the lowest and highest RAFs are al- 
most the same as those of Ceo- The 10 RAFs of C$o cal- 
culated with the B3LYP hybrid DFT method and large 
basis set are in excellent agreement with experiment. Our 
study of Cgo reveals the importance of electron correla- 
tions and basis sets in the ab initio calculations, and our 
calculated Raman results of Ceo demonstrate the desir- 
able efficiency and accuracy of the B3LYP hybrid DFT 
method for obtaining important quantitative information 
on the vibrational properties of these fullerenes. Finally, 
we hope that these calculations will provide an incentive 
for Raman measurements on C48N12. 
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